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In the present paper we consider the BTW model perturbed by random-direction anisotropy with 
strength factor e ranging from to 1 corresponding to BTW and Manna model receptively and 
investigate the properties of the statistical observables for various rates of anisotropy. By increasing 
the e, we observe a cross-over taking place between these models. For small length scales, the curves 
show properties similar to the BTW model whereas in the Infra red limit the corresponding is 
nearly the same as the Manna model. The observations confirm that this perturbation is relevant 
for the BTW fixed point and the infra red limit of the perturbed model is described by the Manna 
model. We also propose a differential equation whose solution properly fits with the the Green's 
function obtained by the simulation. This can help us to obtain the action of the perturbed model. 



I. INTRODUCTION 



In recent years, the concept of self-organized critical- 
ity (SOC) proposed by Bak et. al.*^ has attracted a lot 
of attention as a possible general framework for expla- 
nation of the occurrence of robust power laws in nature 
as it does not require fine tuning of any parameter to set 
criticality. Sandpile models was the first example of these 
systems. The abelian structure of the sandpile model was 
first discovered by D. Dhar and named as Abelian Sand- 
pile Model (ASM)^. Numerous works have been done 
on this model. The connection of this model with span- 
ning trees^, ghost models^, q-state Potts model, Loop 
Erased Random Walk (LERW)^ is known. The differ- 
ent height and cluster probabilities of this model^ and 
its various geometrical exponents are also calculated nu- 
merically and analytically. For a review sed^. Several 
variations of the ASM have been studied in the past with 
a view to understand the parameters that determine the 
different universality classes of self-organized critical be- 
havior. These include models in which particle transfer 
is directed, or models in which the toppling condition or 
the number of sand grains transferred depends on the 
local slope rather than local height^ In this respect, 
it has been realized that stochasticity in toppling rules 
can lead to different critical behavior than models with 
deterministic toppling rules. One of the most interesting 
BTW-variations is the Manna modeP in which after a 
toppling the grains redistribute randomly in a preferred 
direction which is randomly chosen (without dissipation) . 
This model corresponds to a two state model in which 
there is a hard core interaction between two particles in 
the same site that prevents the site to be doubly occu- 
pied. The interesting question then would be what is 
the universality classes which these models belong to. In 
fact, the precise identification of universality classes in 
sandpile models is an unresolved issue. It is generally as- 
sumed that the avalanche size and duration distributions 
follow simple power laws in the infinite-size limit, and 
the departures from such power laws reflect finite-size ef- 
fects. Such effects complicate the estimation of critical 
exponents, since the estimates are sensitive to the choice 
of fitting interval. In the main Manna's papei!^ it was 



claimed that the Manna model lies in the same univer- 
sality class as the BTW model. Real-space renormaliza- 
tion group calculations^^ suggest that different sandpile 
models, such as the BTW and the Manna models, all be- 
long to the same universality class. This result is also 
confirmed by a proposed field theory approach^^ that 
states that all sand pile models are described by the 
same effective field theory at the coarse grained level. 
Using the new exponents, introduced by K. Christensen 
et. al. for the sand pile model^^, A. Ben-Hur et. al. 
showed that these exponents for Manna model are dif- 
ferent from BTW counterparts^^. They claimed that 
the Manna model is in the universality class of random 
neighbor models which is distinct from the BTW uni- 
versality class. Some more exact numerical results also 
confirmed this hypothesis^^. Based on some numerical 
analysis (bias removing), A. Chessa et. al. argued that 
the results of Ben-Hur et. al. are same for the two models 
which imply that they may belong to the same universal- 
ity class^^. They used the finite size scaling (FSS) argu- 
ments in calculating the exponents of size, time and area 
distributions of avalanches of BTW and Manna models 
which is believed that are not single fractals and do not 
fulfil the FSS^^. Many papers argued this result and re- 
ported some exact results about the scaling behavior of 
this modePHH. 

In spite of very much works done on these variations, 
very low attention is paid on the question how these 
models are linked and transform perturbatively to one 
another. Knowing the structure of this transition, one 
can obtain some valuable informations of these models. 
To this end, one can add perturbatively stochasticity to 
the BTW model. It can be done by adding anisotropy 
with random preferable direction with arbitrary strength 
to the toppling rule of BTW model and see how observ- 
ables evolve from BTW to the Manna model. A study 
on patterned and disordered continuous ASM has been 
done in^^ in which it has been shown that the quenched 
disorder lead to an irrelevant perturbation in the confor- 
mal field theory corresponding to the BTW. The plan 
of the present paper is to numerically see how things 
change when the disorder is not quenched. For this, we 
consider the geometrical objects (interfaces) of the per- 
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turbed model and some other statistical functions. The 
properties of the interfaces are directly related to the 
correlation length of the system. In the discrete set up, 
in addition to the correlation length, the system has one 
more characteristic length scale, i. e. the lattice constant. 
In the scales much smaller than the system correlation 
length, and larger than the lattice constant, one expects 
that the statistical features of the curves are properly 
fitted to the unperturbed one. There is an idea to de- 
scribe the geometrical interfaces of 2D critical statistical 
models via growth processes named as SLEP^. It is a 
powerful tool to study the macroscopic interfaces of two 
dimensional systems instead of local fields as is common 
in ordinary CFT. The result of this idea is a complete 
classification of probability measures on random curves 
in (simply connected) domains of the complex plane sat- 
isfying two axioms: conformal invariance and the domain 
Markov propert}^^-^. From the correspondence of ASM 
with the ghost modeP, one expects that this model is 
a c=-2 conformal field theory and with the knowledge of 
the connection between conformal field theory and SLEP^ 
one finds that the ^SM is related to SLE with = 2. 
This test is done in^^ in which numerically is shown that 
the boundary of avalanches of ASM is SLE(2). To gener- 
ate curves running from origin to the infinity from loops, 
one has to cut the loops horizontally and then send the 
end point of the curve to the infinity. In this paper we 
consider perturbed ASM on a cubic lattice and numer- 
ically analyze the properties of interfaces of this model. 
In sections II and III we briefiy introduce ASM and SLE 
respectively. Sections IV and VI also contain some nu- 
merical results from loop statistics and SLE on these in- 
terfaces respectively. 



II. INTRODUCTION TO THE ABELIAN 
SANDPILE MODEL 

Consider the ASM on a two dimensional square lattice 
L X L. For each site we consider the height variable hi 
taking its values from the set {1, 2, 3, 4} which shows the 
number of grains in this site. So each configuration of the 
sandpile is given by the set {hi}. The dynamics of this 
model is as follows; in each step, a grain is added to a 
random site i i.e. hi ^ hi -\- 1^ then if the resulting height 
becomes more than 4, the site toppels and loses 4 sands, 
each of which is transferred to one of four neighbours of 
the original site. As a result, the neighboring sites may 
become unstable and topple and a chain of topplings may 
happen in the system. In the boundary sites, the toppel- 
ing causes to one or two sands to leave the system. This 
process continue until the system reaches to a stable con- 
figuration. Now another random site is selected and the 
sand is released on this site and the process continue. 
The movement on the space of stable configuration lead 
the system to fall in a subset of sets of configurations 
after a finite steps, named as the recurrent states. It has 
been shown that the total number of recurrent states is 



detA where A is the discrete Laplacian. For details se^. 
This model can be generalized to other lattice geometries 
and to off critical set up. For a lattice with d neighboring 
sites, the toppling occurs when hi > then the original 
site will lose d grains and the height of each of its neigh- 
bors will increase by 1. It has been shown that the action 
corresponding to this model is: 



S = / d^z{dOdO) (1) 



where and are complex Grassmann variables. 

Waves; As is mentioned above, the topplings in 
the ASM can be done in any order. One very useful 
way to relax is by a succession of waves of topplings. 
Let the site where the grain is added be O. If after 
addition, O is still stable, the relaxation process is over. 
If it is unstable, we relax it as follows: topple O once, 
and then allow the avalanche to proceed by relaxing any 
unstable sites, without however toppling O again. This 
constitutes the first wave of toppling. If at the end, site 
O is still unstable, we allow it to topple once more, and 
let the other sites relax, until all sites other than O are 
stable. This is the second wave of toppling. Repeat as 
needed. Eventually, site O is no longer unstable at the 
end of a wave, and the relaxation process stops. It is 
easy to see that in any wave, the set of toppled sites 
forms a connected cluster with no voids (untoppled sites 
fully surrounded by toppled sites), and no site topples 
more than once in one wave. (This would not be true if 
the graph had greedy sites). 

Random Anisotropic ASM; One can add ran- 
dom anisotropy to ASM in the following sense; Choose 
random number r = ±1. When the height of a site 
becomes more than 4n where n is some integer, then 
4n grains is transferred to the neighboring sites; n — r 
grains transfer to the 'up' site, n — r grains to the 'down' 
site, n + r grains to the 'left' site and n + r grains to 
the 'right' site. In this respect we define e = ^. e = 1 
and e = correspond to the Manna and BTW models 
respectively. When e = or e = 1 the observables 
have robust power low behaviors up to a characteristic 
length (named as correlation length ^) above which 
the correlation functions falls off rapidly. In a critical 
model, this length is of order of the lattice size and 
when the system size goes to infinity, it diverges. The 
correlation length can be defined as the loop linear size 
i^cut) above which the logarithm of the distribution 
function of gyration radius falls off more rapidly than 
linear. When < e < 1 it is seen that there is another 
characteristic length ^2 in which the behavior of the 
statistical functions smoothly changes. 
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length named as r^^^^, it is seen that N{r) 



where 



(2) 

Tr — 1±0.05. For the lengths higher than r),^^ this distri- 
bution function falls off rapidly, r^^^ may be interpreted 
as the correlation length of the model and in the critical 
case, it is of order of the lattice size. By increasing e, 
another length scale is observed at which the behavior is 
changed. We name this length r^c^t- shown in this 

figure, by increasing e, this point is changed and pushed 
along the origin. This point contains some interesting in- 
formations. Going from small lengths to large scales, we 
see two important lengths. In the small scales the curve 
is locally like BTW's up to the point r^^^^. This suggests 
that the ultra violet (UV) properties of the perturbed 
model is given by the BTW model. In the vicinity of 
r^^^^, the infra red (large scale) limit is reached and the be- 
havior of the graph (more exactly, its slope) is smoothly 
changed to the new one. This and other figures (to be 
shown later) show that in the new regime the properties 
of the model is best fitted to the Manna model (e = 1 in 
the figure). So it seems that the infra red (IR) limit of 
these curves is given by the Manna model and a cross- 
over takes place in between. The slopes of the graphs in 
the IR region is slightly different from the Manna model 
due to the finite size effects. In fact by enlarging the size 
of the lattice, we saw that the slopes in this region get 
closer to the Manna model (it is not shown here). The 
inner graph of this figure shows i.e. the slope of the 
graph for the small scales with respect to e. The hori- 
zontal axis is in logarithmic scale. As is seen, a smooth 
change of behavior takes place when we go from BTW to 
the Manna model. In the vicinity of the e = and e = 1 
the dominant behavior is the BTW's and the Manna's 
respectively and a jump takes place in between showing 
the mentioned cross-over. The same feature is seen in 
FIG [3] in which the distribution function of loop length 
'A/'(/y is shown versus loop length T. In this case, as the 
6 gets non zero values, the deviation from the power law 
(governed on the first part of the graph) takes place at a 
characteristic length depending on e value and this length 
decreases as e increases. For e = 1 the unique power law 
behavior is retrieved N{1) ~ with r^^^ = 1.8^0.1. 

For < e < 1, the IR and UV properties the curves look 
like e = 1 and e = cases respectively. A smooth tran- 
sition for rf from BTW {rj^ = 1) to Manna {r^ = 1.8) is 
seen in the inner graph of this figure. A similar calcula- 
tions was done for the distributions of loop masses and 
the number of waves in the avalanches and loop sizes. 
We observed the same results as above. 
By enlarging the size of the lattice, we saw that this graph 
was totally shifted to the left showing that enlarging the 
lattice size infinitely, the dominant behavior would be the 
Manna's. This tells us that this perturbation is relevant 
for the BTW fixed point and irrelevant for the Manna 
model. For the scales much smaller than the correlation 
length, one expects that the system have a well defined 
fractal dimension. So we computed the fractal dimension 
of waves defined as / ~ r^f and found that for finite e. 



FIG. 1: The schematic view of square lattice. The toppling 
rule is schematically shown (r is the random number taking 
values +1 and -1). 
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FIG. 2: The behaviour of log(N(r)) versus gyration radius 
Log(r), for various dissipations for waves. The inner graph 
shows scaling of their slopes in terms of mass. 



III. NUMERICAL RESULTS; STATISTICS OF 
PERTURBED ASM 



In this section we numerically study the statistics of 
waves and avalanches of ASM to test its dependence on 
e. For this, we have simulated over 100000 independent 
samples and obtained the domain walls between the top- 
pled and untoppled sites on the honeycomb lattice of 
size 2048 x 2048. For the simulation, we considered the 
square lattice. The lattice and toppling rule hase been 
shown in FIG[i]. Consider the wave frontiers of toppled 
sites of ASM. We first study the statistics of the gyra- 
tion radius of loops. In Fig [2] the Log-Log plot of the 
distribution of gyration radius N(r) versus gyration ra- 
dius r is sketched. For the BTW model (e = 0), up to a 
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TABLE I: The r^^^^ and the exponents of gyration radius of 
waves N{r) ^ r^'^, number of waves in an avalanche N{Nw) 
Nw^"^ , the fractal dimension of waves {Log(l)) — Df{Log{r)) 
in terms of e for L = 2048 



IV. DETERMINATION OF THE PERTURBING 
FIELD 



FIG. 3: Log-Log plot of loop length distributions N{1) versus 
/ of waves. The inner graph shows scaling of their slope in 
terms of mass. 



To determine the weight of the perturbing field we 
use two methods. First we directly try to determine this 
weight by using finite size effects. The second method 
we use is the determination of the Green function which 
depends directly to e. 
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FIG. 4: Fractal dimension of wave loops. The inner graph 
shows scaling of this in terms of mass. 



there is a slight deviation from the critical fractal dimen- 
sion (Note that surely for the lengths comparable with 
or larger than the correlation length, fractal dimension 
does not make sense, but to see the behavior of models 
with e's near the Manna model, we calculated the fractal 
dimension for all e's). The interesting feature, as is seen 
in FIG [4], is that this quantity does scale with the log- 
arithm of e in the cross-over region i.e. D^-^^^^ — ~ 
Log(e) and for large e's it does not change and is nearly 
constant near the Manna model. This shows that this 
perturbation is irrelevant for the Manna fixed point. Ta- 
ble [I] shows some informations about the various expo- 
nents of the waves for various e. 



Green function method; In this section using 
Green function of the perturbed model we present a 
method to compute the conformal weight of the per- 
turbing field. As proved by Dhai^ the Green function of 
ASM is defined as follows: suppose that site i is toppled 
by adding a grain. The Green function G{\i — j\) is 
the number of topplings that occur in site j (up to a 
normalization factor) and is proved that the form of this 
function in 2D is G{r) ~ Log(r) which Mathematically is 
the answer of the equation {^dr[rdr])G{r^r) = S{r — r). 
When the action of a critical model is perturbed, the 
action of the model is modified. This modification may 
be of the form: 



n ^ 
k=l •' 



(2) 



in which S** is the action of the conformal field theory cor- 
responding to the critical model (in this case the Eq[l]) 
and S is the perturbed action and \k is the coupling con- 
stants of the perturbations. Then one expects that the 
differential equation of resulting Green function is also 
modified. In our case we can first search to find which 
equation the Green function does satisfy. Then we can 
observe that what is the conformal weight of the per- 
turbing field. In Figjs] the results of the simulations are 
sketched. In this figure the horizontal axis is in logarith- 
mic scale. We see that for low e's the resulting graph is 
logarithmic (up to a characteristic point Rcut which is 
the finite size effect) as expected. But when e increases, 
the graphs do not show this behavior. In the inner graph 
we have shown the Log-Log graphs and see that for e = 1 
the graph is power law G{r) ~ r~^^ with 1. There 
is another way i.e. scaling argument, to study the Green 
function in IR regime. It is known that for a sandpile 
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TABLE II: The explicit amounts of ^(e) for different ( 



model the relation between the exponent of Green f 
tion and the exponent of the distribution function of ^ 
gyration radius is^ 

-Tr^l^df -D = -XI 

where D is the dimension of the space (here is 2) and 

,(2) 

(s) ^ r f where s is the area of the loop and r is its gy- 

(2) 

ration radius. Here =2. In the case e = (r^ 1) it 
is obvious that xi = and the answer will be logarithmic. 
For the case e = 1 also we have (we have calculated above 
that r^^^ 2) xi = 1 which confirms the calculation and 
the simulation. 

In the intermediate values of e we see two different be- 
haviors reflecting their IR and UV properties. For the 
lengths much smaller than the correlation length we see 
that they show similar features like the BTW model 
as expected. For the large scales (lengths much larger 
than the correlation length) they behave like the Manna 
model. This confirms our claim that the IR properties of 
the perturbed models best fit with the Manna model. 
Now we will search for the best differential equation 
which yield the mentioned properties. For large scales, 
the answers of this equation should be like the Manna 
model's i.e. G{r) ~ ^-i±o.i small lengths, 

should have logarithmic form. We have observed that 
the best fit can be achieved by the following differential 
equation: 

^dr[rdrG{r)] - g{e)^dr[rG{r)] = S{r) (4) 

In this equation, ^(e) is some function of e. The numerical 
values of this function has been presented in Table [II]. For 
small values of e, we have ^(e) ~ g0.87±o.04 ^-^^1 use 
this form in the next subsection. 

The Eq.[4] properly yields the IR and UV properties 
of the Green function and for the intermediate lengths 
also fit with the simulation. For example in FIG [6] we 
have shown the simulation and calculated (from Eq[4] 
with g{e = 0.1) = 0.25) results which have been fitted 
well. From the extra term ^dr[rG{r)] we are led to the 
very important conclusion that the conformal weight of 
the perturbing field is 2 A = 1. This statement will be 
further tasted in the next subsection in which we try to 
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FIG. 5: Green function vs Log(r). Inner graph shows the 
dependence of the slope of the first part of curves to the dis- 
sipation. Note that some curves are added by a constant to 
have all the same first point. 
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FIG. 6: Green function vs Log(r). Inner graph shows the 
dependence of the slope of the first part of curves to the dis- 
sipation. Note that some curves are added by a constant to 
have all the same first point. 



yield A from RG arguments. 

RG equation; In the previous subsection we obtained 
the scaling of ^(e) versus e. In this subsection, using 
this construction and the RG arguments, we check the 
validity of this claim. We use the equation that govern 
the off-critical conformal field theory (in vicinity of the 
fixed point). Suppose that the action of the theory is: 
It can be easily proved that the RG equation for this 
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Critical behaviour of the two dimensional statistical 
models can be described by their geometrical features. 
In fact instead of studying the local observables, we can 
focus on the interfaces of two dimensional models. These 
domainwalls are some non-intersecting curves which di- 
rectly reflect the status of the system in question and sup- 
posed to have two properties: conformal invariance and 
the domain Markov propert}/^^. Schramm- Loewner Evo- 
lution is the candidate to analyze these random curves 
by classifying them to the one-parameter classes (SLE^). 
Let us denote the upper half plane by H and 7t as 
the SLE trace i.e. = {z ^ H : Tz < t} and the 
hull Kt = {z e H : Tz <t}. SLE^ is a growth pro- 
cesses defined via conformal maps which are solutions 
of Loewner's equation: 



FIG. 7: Determination of the weight of perturbing field. The 
curve is best fitted to e = 1/1500 



couplings 

dgk 
db 



2(1 - Ak)gk - ^ Cij^kQigj 



(5) 



where Xk = gkO?^^~^^\ is the conformal weight of the 
perturbing field (pk^ b is the length scale, a is the lattice 
constant and Cij^k is the fusion coefficient of the fields 
(fi and cpj i.e. (pi^pi = Xl/c Cij^k^k- In our case, we have 
single coupling constant ^(e) and we can let g{e)'^ 0. 
So we have(^(e) e"^ with n = 0.87 zb 0.04 for smah e's): 



Sg{e)/g{e) Se/e 



n- 



2(1 -A) 



(6) 



where A is the weight of the perturbing field. For us- 
ing the KG equation [RG], we can set a ^ a(l + Sb) or 
equivalently let the size of lattice L = L. So the 

correlation length of the system ^ . To use Eq. jo] , 

it is sufficient to repeat the simulation for L and see for 
which e two cases do match. In Fig[7] we have presented 
the result for L = 2048 and Sb = 0.5 and 62 = 1/3000. 
We see that the best fit is done for ei = 1/1500. Then 
we have 



2A 



Se/e 



0.87^ 
053 



-0.65 0.94 



(7) 



Where e is the averaged value of the e in this interval. 
We have done this test on many such samples with 
various lengths and e's and found the best values for A 
are 0.9 < 2 A < 1.05 in agreement with the obtained 
result in the previous subsection. 



dtgt{z) 



gt{z) - it 



(8) 



Where the initial condition is gt{z) = z and = \^Bt 
is a real valued smooth function. For fixed gt{z) is 
well-defined up to time for which gt{z) = ^t. The 
complement Ht := H\Kt is simply connected and 
gt{z) is the unique conformal mapping Ht ^ H with 
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0{\) as z 00 that is known as hy- 



drodynamical normalization. One can retrieve the SLE 
trace by 7^ = \im^_io g^^ {^t + ^e). There are phases for 
these curves, 2 < < 4 the trace is non-self-intersecting 
and it does not hit the real axis; kt = Jt- This is called 
"dilute phase". But for 4 < < 8, the trace touches it 
self and the real axis so that a typical point is surely 
swallowed as t ^ 00 and Kt ^ jt- This phase is called 
"dense phase". However, there is an important property: 
The frontier of Kt i.e. the boundary of Ht minus any 
portions of the real axis is a simple curve and is locally 
SLE^ with k = This duality links models in dilute 
phase to one model in the dense phase and vice versa 
e.g. the ASM {k, = 2) to the Uniform Spanning Tree 
(UST) {k: = 8). The main question "what is the relation 
between SLE and CFT" is answered by M. Bauer 
and D. Bernard^^. They showed that the boundary 
condition changing (bcc) operator in SLE correspond 
to a degenerate field with a vanishing descendant at 
level two and conformal weight hi- 2 

central charge = '^^~^2^f • 



^ in CFT with 



SLE Out of criticality; Now consider the system 
out of criticality. In this case the conformal invariance 
of the system is broken and the the system correlation 
length C will have a crucial role in statistical properties 
of the random curves. So if we apply the Loewner 
uniformizing map, the resulting domains are not more 
equivalent due to absence of conformal invariance. At 
scales much smaller than the correlation length, i.e. in 
the ultraviolet regime, the deviation from criticality 
is small, and the interface should look locally like the 
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critical interface. This means that over short time 
periods, the off-critical driving function should not 
be much different from its critical counterpart. In the 
other hand, at large scales (with respect to ("), i.e. in 
the infrared regime, the interface may look like another 
SLE with a new hZir. The reason is; when we integrate 
out the small distances to reach the large distance 
properties, the regions which is formed by SLE trace 
for lengths smaller than the correlation length (with 
diffusivity /^t/y), may be seen as points that the SLE 
trace with the new diffusivity (hCir) corresponding to the 
new infra red fixed point crosses. The example is Ising 
model. At criticality = 3^ but if the temperature is 
raised above the critical point, renormalization group 
arguments indicate that at large scale the interface looks 
like the interface at infinite temperature i.e. percolation 
with Kir = d^. 

Now consider a curve that starts from origin and end 
on a point on real axis (xqo)- Then by using the map 
^ = ^00^/(^00 — ^)7 one can send the end point of 
the curve to the infinity. In this respect, the function 
ht = (j)ogtO(j)~^ describes chordal SLE. It is easy to show 
that the equation governing on ht is: 



to Xoo and by demanding that this semicircle involve Zi 
one obtain: 



dtgt = 2/{{4>{gt){4>{9t) - it)])- 



(9) 



This mapping is not hydrodynamically normalized i.e. it 
does not fix the infinity, instead it fixes the ending point 
of the curve. 



VI. NUMERICAL RESULTS; 
SCHRAMM-LOEWNER EVOLUTION 

In this section we present some numerical results ob- 
tained by applying SLE on the critical and off-critical 
abelian sindpile model. The the frontier of avalanches 
form the set of loops with discrete points. As in the 
chordal SLE, the curve goes from a point of the real axis 
to a point in the infinity and here we have loops, we are 
to use a trick to generate desired curves. Having these 
loops, one can cut them with a straight line to generate 
interface curves starting from the origin and ending at 
some point on the real axis (xoo). Then using the map 
(j) = Xooz/{xoo — z), that fixes the origin and sends the 
ending point (xoo) to the infinity, we will have a curve on 
the upper half plane. Then by applying the chordal SLE 
formalism and the proper uniformizing map step by step, 
one can emerge for these discrete curves. The essen- 
tial assumption is that is partially (in each interval) 
constant, then it can be easily proved that the mapping 
that can be used to uniformize the curves is^^: 



Gt{z) 
where r] 



{xoo - z) +a(^^,z) 
^-^(0 and a{it,z) 



(10) 



Re{zi)Xoo - [Re{zi)f - [Im{zi)f 
Xoo - Re{zi) 



t = 



Re{z,)'x^ 



4{Re{zi) - Xoo)^ ^ Im{zi)Y 



(11) 



(12) 



^/x^{z - r]Y + 4t(xoo - ^)^(^oo - v)- This map- 
ping uniformizes a semicircle that is extended from r] 



FIGjs] contains the graph (^|) — (^t)^ versus t for the 
critical case. As it is explicit in the graph, has the 
expected behaviour: (^t) and ((^|) — (^t)^) = nt 
with n — 2.0 ± 0.1. We note that the initial portion 
(0 < t < 1000) of the graph is different from the 
remaining and has been ignored. The reason is that for 
these times, the effect of finite size (lattice constant) on 
the curve growth is important as the size of the curve is 
comparable with it. So in the remaining of the paper we 
will ignore this portion. We however are not concerned 
about the effect of the size of the system, because it take 
very long time for such a fractal to have a linear size of 
the system order. 

For the off-critical model, as stated in Sec. [V| we have 
two important scale limits. For small scales (scales 
much smaller than r^^^^) the interface should look locally 
like the critical interface at UV fixed point. At large 
scales however, the interface may behave like a SLE 
corresponding to the IR fixed point with Kir . We 
consider the perturbed ASM as described bove and 
analyze the resulting driving function. The important 
quantity which can be extracted from driving function is 
K which in the critical case is obtained from the relation 
(^?) ~ ^^t)^ — III off-critical case we may observe 
two slopes for the graph: one for UV region (the result- 
ing K should not be much different from the critical one) 
and another for IR region. In between the curve may 
have complex behaviors. In FIG [9] (^^) — (^t)^ versus 
t is shown for some e's. The slope of the graph for the 
critical (e = 0) one is k = 1.95 ± 0.1 in agreement with 
other numerical results^^. When e becomes non-zero, 
the graphs does not show simple linear behavior. As is 
seen from this figure, there are two transition points: 
The first one is the earliest time in which the graph 
separates from the critical one i. e. is the first transition 
point from UV to the "cross over" region (in the FIG [9] 
is shown as Ti). The next (T2) is the transition from 
the "cross over region" to the the IR region. In this 
region, the slope is the same as the Manna's i. e. 
e = 1. These transition points depend on e and increase 
as e decreases and in the case e = become infinite. We 
can investigate the behaviour of the curves at UV and 
IR regions well separated from the crossover region. The 
result is that the curves are linear in each region with the 
same slope in each region. These feature for IR regime 
has been shown in Fig [9] and magnified in FIG 10 in 
which is seen that all have the same slope Kir — 1.65 ±0.1. 
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FIG. 8: The averaged squared versus t shows the diffusivity 
= 2.0 ±0.1. 



n 8=1 

o c = 1/1 00 

A e= 1/500 
e = 1/700 
• e= 1/1 000 



K.=1.65 



FIG. 10: Kir for various dissipations in region II. 
VII. CONCLUSION 
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FIG. 9: {{^t) — i^t)^) versus t for different dissipations. 



In this paper, we analyzed the statistics of wave and 
avalanche frontiers of continuous random anisotropic 
ASM. The BTW model corresponds to the perturbation 
parameter e = and Manna model to e = 1. It has 
been shown that a cross over takes place between these 
two models. We studied the behavior of some statisti- 
cal observables and found the conformal weight of the 
perturbing field by two methods: Green function and 
RG arguments. Each of them confirm that the weight 
X = 1. Using SLE for the geometric curves of the per- 
turbed model we showed that there are two important 
length scales in which the corresponding SLE parameter 
'/t:' is different. These scales are determined with respect 
to the correlation length. Using SLE, we found numer- 
ically that for the scales much smaller than the correla- 
tion length, the curves have the same properties as the 
UV critical model (BTW) with nearly the same K:. For 
scales much larger than it, also we found that the curves 
acquire the new kir — 1.65. 



^ P. Bak, C. Tang, and K. Wiesenfeld, Phys. Rev. Lett. 59, 
381 (1987). 

2 D. Dhar, Phys. Rev. Lett. 64, 1613 (1990). 
^ S. N. Majumdar and D. Dhar, Physica A 185 (1992) 129. 
^ S. Mahieu, P. Ruelle, Phys. Rev. E 64 (2001) 066130. 
^ S. Majumdar, Phys. Rev. Lett. 68, 2319-2331 (1992). 
^ S. N. Majumdar and D. Dhar, J. Phys. A: Math. Gen. 24 
(1991) L357. 

D. Dhar, Physica A 369 (2006) 2970. 
^ S.S. Manna, J. Phys. A 24 (1991) L363. 
^ L.P. Kadanoff, S.R. Nagel, L. Wu, S. Zhou, Phys. Rev. A 

39 (1989) 6524. 

S.S. Manna, Physica A 179 (1991) 249. 

Y.-C. Zhang, Phys. Rev. Lett. 63 P. 470, (1989) 

D. Dhar, R. Ramaswamy, Phys. Rev. Lett. 63 (1989) 1659. 

A. Ben-Hur, O. Biham, Phys. Rev. E53, N2, R1317 (1996). 

K. Christensen and Z. Olami, Phys. Rev. E 48, 3361 



(1993). 

A. Chessa, H. E. Stanley, A Vespignani, S. Zapperi, Phys. 
Rev. E 59, Nl, R12 (1999). 

S. Lubeck, arXiv:[cond-mat.stat-mech] /9910374vl (1999). 
^'^ R. Dickman, J.M. M. Campelo, arxiv:cond-mat/0301054vl 
(2003). 

H. N. Huynh, G. Pruessner, and L. Y. Chew, arxivxond- 
mat/1106.0406vl (2011). 

L. Pietronero, A. Vespignani, and S. Zapperi, Phys. Rev. 

Lett. 72, 1690 (1994); A. Vespignani, S. Zapperi, and L. 

Pietronero, Phys. Rev. E 51, 1711 (1995). 

R. Dickman, A. Vespignani and S. Zapperi, Phys. Rev. E 

57, 5095 (1998); A. Vespignani, R. Dickman, M. A. Munoz, 

and S. Zapperi, e-print cond-mat/9806249. 

N. Azimi-Tafreshi, S. Moghimi-Araghi, arXiv:0906.0459vl 

[cond-mat.stat-mech] 2 Jun 2009. 

H. Saluer, B. Duplantier, Phys. Rev. Lett. 58. (1987) 2325, 



9 



A. Coniglio, Phys. Rev. Lett. 62 (1989) 3054. 
Oded Schramm: Scaling limits of loop-erased random 
walks and uniform spanning trees. Israel J. Math., 
118:221288, 2000. 

M. Bauer, D. Bernard. Conformal field theories of stochas- 
tic Loewner Evolutions. Comm. Math. Phys., 239(3) ;493 
521, 2003. 

A. A. Saberi, S. Moghimi-Araghi, H.Dashti-Naserabadi, S. 
Rouhani, arxiv:0812.0939vl, cond-matt.stat-mech. 
D. Bernard, G. Boffetta, A. Celani, and G. Falkovich, PRL 
98, 024501 (2007). 



Michel Bauer, Denis Bernard, Kalle Kytola, 
arxiv:0712.1952vl,math-ph(2007). 

J. Cardy, arXiv: 0503313v2 cond-mat.stat-mech (2005). 
J. Cardy, "SLE for Theoretical Physicists", arXiv:cond- 
mat/0503313v2 (2005). 

D. V. Ktitarev, S. Lubeck, P. Grassberger, V. B. Priezzhev, 
Phys. Rev. E 61, P. 81 (2000). 

G. Mussardo, 'Statistical Field Theory', Oxford University 
Press, (2010), sec. 15. 



